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Abstract 

We analyze a class of dynamical systems of the type 

d n (i) = c n _ia n _i(t) - c n a n+1 (t) + /„(*), n£N, a = 0, 

where /„(i) is a forcing term with f n (t) 7^ only for n < n* < 00 and the 
coupling coefficients c n satisfy a condition ensuring the formal conservation of en- 
ergy I J2 n a n (t)| 2 . Despite being formally conservative, we show that these dynami- 
cal systems support dissipative solutions (suitably defined) and, as a result, may admit 
unique (statistical) steady states when the forcing term f n (t) is nonzero. This claim is 
demonstrated via the complete characterization of the solutions of the system above for 
specific choices of the coupling coefficients c n . The mechanism of anomalous dissi- 
pations is shown to arise via a cascade of the energy towards the modes with higher n; 
this is responsible for solutions with interesting energy spectra, namely E|a n | 2 scales 
as n~ a as n — > 00. Here the exponents a depend on the coupling coefficients c n and 
E denotes expectation with respect to the equilibrium measure. This is reminiscent of 
the conjectured properties of the solutions of the Navier-Stokes equations in the in- 
viscid limit and their accepted relationship with fully developed turbulence. Hence, 
these simple models illustrate some of the heuristic ideas that have been advanced to 
characterize turbulence, similar in that respect to the random passive scalar or random 
Burgers equation, but even simpler and fully solvable. 



1 Introduction and main results: Life starts after blow-up 

So little is understood about hydrodynamic turbulence that there is not even consensus 
on what it is. However, most physicists would agree on the following heuristic picture 
which has emerged from the works Kolmogorov, Onsager, Richardson, etc |Fri95 1. In 
this picture, (fully developed) turbulence refers to the idealized state of an incompress- 
ible fluid described by the Navier-Stokes equations in the limit of vanishing molecular 
viscosity. In this limit, the Navier-Stokes equations formally reduces to the Euler equa- 
tions, and the turbulent solutions should be the most regular solutions of the Euler 
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equations which dissipate energy. This is referred to as anomalous dissipation and is 
best visualized in the Fourier representation. There it corresponds to a cascade of en- 
ergy from the small wavenumbers (large spatial scales) where energy is injected (either 
via the initial condition or by a forcing term in the equation) towards larger and larger 
wavenumbers (smaller and smaller scales), up to infinity where energy should eventu- 
ally be dissipated. It is also believed that the cascade of energy implies that the energy 
spectrum of the turbulent solutions have a power law decay in the wavenumber whose 
rate can be deduced by dimensional analysis and is | in three-dimension of space. 

Turbulence theory (as we shall refer to the heuristic picture above) also discusses 
more advanced and more controversial topics such as intermittency. But, without even 
going into those, most mathematicians would agree that a rigorous confirmation of the 
basic predictions of turbulence theory is already a tremendous challenge. The best 
known results on the Navier-Stokes and Euler equations which corroborate the above 
were obtained in I C ET94I |EyiO 1 1 [DROP | . These works only indicate that turbulence 
theory is not blatantly inconsistent. Simpler models, such as randomly forced Burgers 
equation or Kraichnan's model of passive scalar advection (see e.g. IE01I IFGV01I 
for reviews), have also been used to demonstrate that parts of turbulence theory make 
sense in terms of anomalous dissipation of the weak solutions of the inviscid Burgers 
equation and the spectrum of energy of the solutions that this implies. Even these 
simple models remain surprisingly complicated to analyze and a full characterization 
of the statistical properties of their solutions is still lacking. 

One of the purposes of the present paper is to illustrate turbulence theory on even 
simpler models. Many (if not most) of the realistic features have been neglected in our 
models. Yet, the models possess a rich range of behaviors which depend on the details 
of the interactions. They provide a simple class of exactly solvable models which can 
be useful in understanding the inner workings of some energy transfer mechanisms. 
The solutions of these models are also consistent with much of the claims of turbulence 
theory. In a way, they offer a setting for the skeptical mathematician to understand the 
motivation behind these claims, and if this paper succeed in doing this, we will have 
achieved our main goal. 

Next, we introduce the models that we will investigate and we summarize the prin- 
cipal results of the paper. As we will see, the most interesting and meaningful solutions 
of these models are solutions which have blown-up, such that they have become infi- 
nite in some norm. This justifies our claim that "life starts after blow-up": Disregarding 
these solutions as nonsensical, as one may be tempted to do at first sight, would, in fact, 
completely miss the most interesting phenomena displayed by the models. 

1.1 A linear shell model 

Consider the equation 

a n (t) = c[(n - I)a„_i(£) - na n+1 (t)] 

for n E N with the boundary condition ao(t) = for all t. If c > 0, we can rescale time 
to fix c = 1; observe also that if a n satisfies the equations with the parameter c < 
then a n (i) — (— l) n+1 a n (t) satisfies the equations with parameter |c|. 

In light of these considerations, we set c = 1 and focus our attention on 

a n (t) = (n - l)a„_i(f) - na n+ i(t) (1.1) 

for ri 6 N with the boundary condition a (t) = for all t. 
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Although we will see that this calculation is not always correct, on the formal level 
one has that 

1 d x 

2^ X! l a «W| 2 = ~ V a n-l(t)a n (t) - na n+1 (t)a n (t)] 

n=l „=1 

OO oo v 7 

= ^2 na n {t)a n+1 (t) - na„(t)a n +i(t) = . 

n—1 n—1 

The second equality is only formal as it assumes that the sum Xm=i Tla n(Qin+i(i) is 
finite and absolutely convergent. To understand this further, consider the evolution of 
the partial sum Y, n <N l a n(0| 2 - For N £ N, 

\j t J2 l a » (i) ! 2 = -Na N (t)a N+1 (t). (1.3) 

n<JV 

The validity of i ll ,2b necessitates 

lim Na N (t)a N+1 (t) = 0. (1.4) 

If this condition is not satisfied, then the formal manipulation in II. 2> does nof hold and 
the seemingly conservative coupling term in II. Il l may become a source of anomalous 
dissipation. We make the concept of anomalous dissipation precise in Section|2] But, 
roughly speaking, it is when seemingly conservative terms have a dissipative effect on 
the system. 

In the context of equation II. 1> . anomalous dissipation seems to require that the 
limit as N — > oo of the right hand side of II. 3i be negative. In other words, equation 
J 1 . 1 1 is dissipative at time t if 

liminf na n (t)a n +i(t) > 0. (1.5) 

n — >oo 

If we make the reasonable assumption that lim„ a n +±/a n £ (0, oo), then from II. 5> 
the solution of equation II. 11 will be dissipative with a finite dissipation rate provided 
that ^ 

a n (t) ~ —= as n — * oo. 1 (1-6) 

Vn 

At this point some readers may be skeptical since one typically considers equations 
like ( ll.lt with initial data in £ 2 , the space of square-summable sequences. However, 
we will see (Theorem B.ll in Section|3j that equation J 1 . 1 1 has solutions which exist for 
all time provided 

limsup|a n (0)| 1/n < 1 . 

n — >oo 

This condition admits a large class of initial conditions including those which scale 
as (II. 6> . In Section|3] we will also see that J 1 . 1 1 possesses a wide verity of behavior 
including conservative, dissipative, and explosive solutions. 

It might be tempting to dismiss these non-conservative solutions as non-physical 
solution arising from pathological data. We now discuss why is not the case. 

'We say that f(m) ~ g(m) as m — > oo if their exists an ni\ and c > 1 so that if m > m,\ then 
ig(m) < /(m) < cg(m). Similarly, we say that f(m) ~ g{m) as in — > mo if lim f(m)/g(m) = 1 as 
m — > mo. 
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Consider equation (ll.lt with a white-noise forcing in the first coordinate: 

ajt) = (n - l)a n _i(t) - na n+1 (t) + l n=1 W(t). (1.7) 

where l n=m is 1 if n = m and otherwise, and W(t) denotes a standard Brownian 
motion, i.e. Gaussian process with mean zero and covariance KW(t)W(s) = min(i, s). 
If one were to accept the formal calculations in il.2\ . showing energy conservation, 
then 

oc oo 

E^|a„(t)| 2 =E^|a„(0)| 2 +i 

n—1 n—1 

if the energy is initially finite. Hence, in the forced system energy seems to grow 
linearly with time and at t = oo one expects the system to have infinite energy. These 
solutions which "blow-up" (in the sense that they have infinite energy) are the most 
interesting and relevant. In light of the discussion above, one might expect that the 
energy of the system would grow to be infinite and arrange the a n so that the calculation 
in II. 2> is not valid since the sum is not rearrangeable. Onsager would then predict 
that the system would evolve to the state in which the a n decayed as fast as possible 
but sill dissipated energy in the sense that (II. 5> holds. The reasoning which leads 
to II. 6> strongly suggests that the \a„\ should scale as 1/^/n. In fact, if the system 
is to reach some equilibrium the effect of the dissipation must exactly balance that 
of the forcing. Specifically, in the stochastic setting when the forcing is W(t)l n —i, 
lim n nE(a n (()fl n -n(i)) — * 1 as t — > oo. 

All of these conclusion turn out to be correct. In particular, in Section[6]we prove 
that if J2 n l a n(0)| 2 < oo then the solutions converge to a unique random variable 
a** = (a**, o", • • • ) which is Gaussian with mean zero and whose distribution is the 
unique stationary measure for the system. Furthermore, this equilibrium state has a 
structure which is consistent with the anticipated \j \fn scaling: 

lim Ea n (t)a m (t) — Ea**a** — . (1-8) 

t— >oo n + to — 1 

We also show that similar behaviors are observed with a different type of forcing. In 
particular, if the forcing is constant, 

d n (t) = (n - I)a n _i(t) - na n+1 (t) + l n=1 , (1.9) 

then the system evolves to a unique steady state a* which scales as 

lim a n (t) = < = ~ J^- . (1.10) 

In summary, we see that the forced systems' energy grows linearly with time if the 
energy is initially finite. Asymptotically, the system rearranges itself so that it reaches 
a state which dissipates energy at t = 00. This state is chosen so that the dissipation 
rate matches the energy flux into the system from the forcing. Since the energy flux is 
finite, this leaves | a n | ~ 1/ %fn as the only choice. A slower decay rate would produce 
an infinite rate of dissipation and a faster decay rate would produce a system which 
conserved energy since the limit in (II .51 would be zero. 
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1.2 A second linear shell model 

We introduce a second model which also exhibits interesting but different "blow-up" 
behavior. Consider 

b n = (n - l)(n - \)b n -i - n(n + \)b n+x (1.11) 

for n 6 N with the boundary condition bo(t) = for all t. As in the previous subsec- 
tion, the unforced equation formally conserves energy since 

1 A 00 00 

2^I^» (t) = S [(n - 1)(n ^ )6 "- l(t)6 '^ (L12) 

n—1 n—1 

This equality (as in the case of equation {Q) ) is only formal since, in general, the sum 
cannot be rearranged. As before, to gain insight we consider the partial sums. For 

N e N, 

\j t Y. b n ( V = + §)M*)&iv+i(*). 

n<JV 

With this in mind, dl.il> will be called conservative if 

lim N(N + \)b N (t)b N+l (t) = 0, 

and dissipative if 

liminf N(N + ±)b N (t)b N+1 (t) > 0. 

A^oo 

Unlike the case of equation Jl.lt . if one assumes that liniAr^oo 6jv+i/^jv exists and 
is in (0, oo), then the solution of equation II. Ill will be dissipative if b n ~ 1/n as 

n — > cx). 

A solution satisfying &„(£) ~ 1/n (if one exists) has finite energy: J^Li &n(*) < 
oo. Thus, this model differs from the first example in that the system can dissipate 
energy even when the total energy is finite. 

While we do not prove a general existence result as broad as for II. H . we do show 
(in Theorem l9.1> that the dynamics for ll.ll> are well defined if 

oo 

^(-1)"M0) < oo. (1.13) 

71=1 

This is sufficient for our needs: In particular, it covers the case when b n ~ 1/n as 

n — > oo. 

The differences between the first and second models are greater than simply the 
scaling. When started with initial conditions having finite energy the first model con- 
serves energy. In fact the regularity at time t is the same as the regularity of the initial 
condition. In contrast, if we start II. 1 li with initial data satisfying I1.13> (and hence 
with finite energy), the energy decays with time. Furthermore, for almost every t > 
one has that b n (t) ~ 1/n as n — > oo and there exists T, depending on the initial data, 
so that if t > s > T then 

oo oo oo 

E < E ^ < E ^ < °° 

n—1 n=X n—1 
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and t>n(t) — * and t — > oo. Turning to the forced setting, consider 



(n — l)(n 




n(n + \)b n+1 + /(i)l n =i. 



When f{t) = Wit) then 




1 



as 



n 



oo. 



n(n + m) 



If fit) = 1> we nave 



h* 

n 



l 



as 



n 



oo. 



n 



1.3 Inviscid limits of the first model 

In practice, one is often interested in understanding the limit of equations when the 
explicit sources of dissipation are removed. To explore this question we investigate 
equation (II. II with the addition of an overtly dissipative term and then study the limit 
as the dissipation is removed. 

To understand our motivation, recall that we have seen that if the first model is 
started with finite energy initial data then the formal calculation presented in II. 2> is 
valid for all finite times as the energy must be infinity for II. 2> to fail. Yet, as time tends 
to infinity, the forced system converges to a steady state with infinite energy for which 
the calculation presented in ( II .21 fails. In contrast, in the second model the analogous 
calculation, given in 11.121 . fails at almost every positive time since bit) ~ 1/n as 
n — > oo for almost every t > 0. 

Since the coupling term produces dissipation at finite times in the second model, 
it is most interesting to study the effect of extra, explicit dissipation in the first model. 
To this end, we consider the stochastically forced version of the first model with extra, 
explicit dissipation sufficient to keep expected energy of the system finite for all times. 
In particular, the calculation in II. 2> is valid in the equilibrium state. We are interested 
in the structure of this invariant state and how it converges to the steady state without 
the explicit dissipation (as the dissipation is removed). 

We will consider two cases, one for which the dissipative term is lower order than 
the coupling term and one for which it is higher order. Specifically, forp £ {0, 1} and 
any v > 0, consider the equation 

a nM = -2vin - Vfot n , v + [in - l)a„_i^ - na n+lM ] + l n=1 Wit). 

The case p = corresponds to the lower order damping and is analogous to what is 
called Eckman damping in the context of fluid mechanics. When p — 1 the perturbation 
is higher order than the coupling term and behaves as a viscous term in the language of 
fluid mechanics. As in the previous examples, we force the first coordinate with white 
noise. Assuming that \ a n,A0)\ 2 < oo, it is straight forward to see that \a UtV it)\ 2 
stays finite and uniformly bounded in time for any v > 0. Hence, the system remains 
conservative for all time. Furthermore, as t — > oo the system converges to a random 
variable a** v whose distribution is the unique stationary measure for the system. Direct 
calculation in the spirit of II. 2> shows that 



oo 



1 



n=l 
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Thus, there is no anomalous dissipation in the system: All of the dissipation which 
balances the forcing comes from the term —2v(n — l) p a n (t). In section[8] we will see 
that for p £ {0, 1}, 

E[<*„ - a**] 2 ^ as v -> 0. 

Both of these steady states are Gaussian with mean zero. The way in which the 
variance of converges (as v — > 0) to that of the v = steady state a** is different 
for the two values of p we consider. When p = one has 



This shows that variances still decay like a power of n. Notice that for all v > the 
total energy is finite. When p = 1 we do not obtain an exact formula but rather that 

<E[<*] 2 < — — — (1.15) 



(« + v) 2n + 2 2n + 1 - L n ' vS ~ K 2n 2n 

where k 2 = 1 + v 2 . Since k > 1 when i/ > 0, the a**„ behave as the limiting a** for 
small n but decay exponentially for large n. 

1.4 Organization 

The remainder of the paper is concerned with proving the statements made in this sec- 
tion. Section|5]contains a precise discussion of the concept of anomalous dissipation. 
In Section [3] we return to the first of the two models introduced in Section [2 and il- 
lustrate the range of possible dynamics by considering specific initial conditions for 
which the system can be explicitly solved. In Sections|4]and|5] we prove the existence 
of solutions and describe their properties for a wide range of initial data. In Section^ 
we give the proofs of all of the preceding results. Section[6]discusses the forced setting 
for the first example and Section|8]discusses its inviscid limit. Finally, in Section[9]we 
discuss the second model introduced in Section^ we first describe the qualitative be- 
havior of solutions; then, we prove existence and uniqueness results both in the forced 
and unforced situations. 

We note that some of the results about anomalous dissipation and the scaling of 
the solutions can also be obtained by formally taking the continuous limit in n of Jl.lt 
and (11.1 II . In this limit, these equation formally reduce to hyperbolic conservation 
laws which were analyzed in |SriQ5]. 



2 Preliminary: Definition of anomalous dissipation 

The concepts of energy conservation, dissipation, and explosion are straightforward 
when the total energy of the system is finite. A system is conservative if the energy 
does not change with time. A system is dissipative (or explosive) if the total energy 
decreases (or increases) with time. 

However, as the example in the previous section showed, it is possible to have 
solutions which one might call dissipative even though the total energy is infinite. We 
give a definition of the above terms which can be applied to situations where the total 
energy is infinite. 



The rich behavior of the first model 



8 



Given a time-dependent sequence, {a n {t)} ne ^, define the energy in the block M 

to N, M < N, by 



JV 



n=M 

A given block £m.n is dissipative at time t if £m,jvO0 < 0. Similarly, we will say it 
is explosive if £M,N(t) > 0. If £M,N(t) = then we say the block is conservative at 
time t. If £M,N(t) — for all M, N, and t, then the system is at a fixed point. (Note 
that this is consistent with Example l3.5l in the next section.) 

We will say that the system is locally dissipative (locally explosive, or locally con- 
servative) at time t if every finite block is dissipative {explosive or conservative) at time 
t. 

In contrast, we will say that a system with £o,oo = oo is dissipative at time t if the 
limit ^ 

limsup — £o.N(t) < 0. 

n^od at 

We say it is explosive at time t if 

liminf —£q N(t) > 0. 
jv^oo at 

When the limit exists we will refer to its absolute value as the rate of energy dissipation 
or the rate of energy explosion depending on the inequality which is satisfied. We say 
that the system with £q,oo = oo is conservative at time t if 

lim ^r-£o,N(t) = 0. 

n^oo at 

If we do not state time explicitly for any of these property, we mean that the properties 
holds for all finite times. 

As the examples of the next section show, it is possible for linijv £q,n = °° while 
liniTv £q,n — c < 0. 

Remark 2.1 It is important to notice that the above categorizations are not exhaustive. 
It is possible for a system not to fit into any of the categories. This is only an issue when 
the energy is infinite as we use the definitions at the start when the energy is finite. 



3 The rich behavior of the first model 

The system given by ( II. U possesses a number of interesting properties beyond those 
listed in the introductory section. In this section we explore the behavior through a 
number of examples. A relatively complete theory of the equation will be given in the 
two sections which follow. We begin with an existence result which covers all of the 
examples presented. 

For an infinite vector a — (at, a%, ■ ■ •) with a, £ R, define p(a) by 

-i- = limsup . (3.1) 

p{a) n~>oc 

(p(a) is simply the radius of convergence of the power series J2 n a « z ™)- The following 
theorem gives an existence and uniqueness result for sufficient for present needs. 
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In particular, it allows initial data with infinite energy |a„(0)| 2 = oo). A more 
complete description will be given in Theorem l4.1l of Section|4] where we describe the 
methodology for solving II. 1> . 

Theorem 3.1 If a(0) = (ai(0), a2(0), • • • ) is an infinite vector of initial conditions such 
that p(a(0)) > 0, then there exists a unique solution a(t) to ll.lt with initial conditions 
a(0) which exists at least up to the time = arctanh(p(a) A 1). 

This existence result, whose proof is given in Section^ covers a wide class of initial 
data. The dynamical behavior of II. Q is quite rich. We list a number of exact solutions 
which display the range of possible behaviors. Explanations of how these results are 
obtained will be given in section l4~Tl A general discussion of the qualitative properties 
of solutions of II. li will be given in Section[5] 

Example 3.2 An energy conserving pulse heading out to infinity: Fixing ai(0) = 1 
and a„(0) = for all n = 2, 3, . . . results in the dynamics 

tanh"- 1 ^) 

Qn ^ = U7T~ ■ 

cosn(i) 

Even though the solution decays to zero pointwise in n as t — > oo it conserves energy: 
Em=i l a "W| 2 = S^Li l a n(0)| 2 = 1 f° r a ^ t > 0. The fact that it conserves energy 
is consistent with the observations in equation II. 2> and ll.4t because 



, V I '1(111 

\a n (t)\ 2 = -Na N (t)a N+1 (t) = — 



Id ^ 2 A7 /_t_\ ,^ Atanh^CQtanri (Q 



2dt- ^™ ~ "" flW " flW " co S h 2 (i) 

n<N v 7 

as A — > oo. The dynamics of this solution can be understood as pulse moving out to 
larger and larger A with time while simultaneously spreading out. A simple calculation 
shows that a,N(t) reaches its maximum at a time asymptotic to \ log(l + 4 A) as A — > 
oo. Hence, as t — ► oo the n for which a n which is cresting at time t scales as 4 exp(2i). 

Example 3.3 Dissipative solution with finite dissipation rate: If 

(2n)! 1 



x„+i(0) = = for n = 0, 1, 

2 z "(n!)^ V7rn 



then 



m e ~ t>2 V- (2(n - m))\(2m)\ 
Vcosh(i) ^— '„ 2 zn ((n — mnmlY 

For each fixed n, we have 

a n (t) ~ \/2e _t as £ — > oo, 

so that the solution decays to zero pointwise in n as t — > oo. On the other hand, since 
for any fixed time t 

e"'/ 2 1 
a n (t) ~ -^^^=^^= asn^ oo, 
Vcosht V 71-71 
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it follows that J2^=i \ a n(t)\ 2 = +°° f° r a ^ ti me t >0 and 

\t + M*>l 2 = -Na N (t)a N+1 (t) -> e —— < 

2 at ' 7r cosh t 

n<N 

as JV — > oo. For this solution, the calculation in II. 2t does not hold and the above 
inequality can be interpreted as a form of anomalous dissipation with finite dissipation 
rate. 

Example 3.4 Dissipative solution with infinite dissipation rate: If a„(0) = 1 for all 
n = 1, 2, then 

a„(i) = e _t for all n = 1, 2, . . . 

This solution decays to zero pointwise in n as t — > oo and X^Li | a n(£)| = +°o for 
all < > 0. Notice that 

~ l«™W| 2 = -^w(*)flA'+i(<) = -Ner 2t -» -oo 

n<AT 

as iV — ► oo. The formal calculation in i ll ,2b does not hold for this solution and, in 
terms of the definitions of Section we view this as a form of anomalous dissipation 
with infinite dissipation rate. 

Example 3.5 A fixed point: If a2 n (0) = for n = 1, 2, . . . and 

(2n)\ 1 

a 2n +i(0) = 2 2 = for n = 0, 1, . . ., 

2 z ™(n!) z .y%n 

then 

«n(Q = fln(0) for all 71 = 1, 2, . . . 
This solution is a fixed point of II. 1> . Notice that 

Id N 

2dt^ |a " (t)|2 = - Na N(t)a N +i(t) = 

n— 1 

for all N S N, but ^Z^Lj l a "(Q| 2 = +°° smce c&2n+i ~ 1/y/nfi as n — > oo. 

Example 3.6 Explosive solution with infinite explosion time: If a„(0) = (— 
for n = 1,2,.... then 

a n C9 = (-l)" +1 e* for all n = 1, 2, . . . 
In this case, Y^™=i l a «W| 2 = +°° an d 



Id N 



2dt 



^ M*)| 2 = -iV aA r(i) ajv+1 (t) = iVe 2 * -> +00 



as TV — > oo. Thus, II. 2t does not hold for this solution and we see that the above limit 
is consistent with an infinite explosion rate for the solution. This is consistent with 
Theorem O.ll as p{a) = 1 and t* = arctanh(l) = oo . 
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Example 3.7 Explosive solution with finite explosion time: If a„(0) = (— l) n+1 a™ 
with a > 1 for n = 1,2,..., then 

{-l) n+1 a n 

= for t < <* = arctanh(l/a!) and all n = 1,2, . . . 

cosn(i) — a sinn(i) 

This solution blows up at t = t*. In this case, X)^Li l a «(^)| 2 — +°° f° r a U t < t* and 

1 rf A r a 2Ar+1 

9^7 H l fl " (i) ! 2 = -Na N (t)a N+ i(t) = — , 2 -> +oo 

J «* ( cosh(t) - a smh(i)) 



as N — > oo for all t < t„. Notice that this example is consistent with Theorem O.ll as 

p(a) = 1/a < 1. 

The above examples demonstrate the rich range of behavior of the model. In partic- 
ular, some solutions grow coordinate-wise in time while others decay. The following 
result gives a criteria for the later. 

Theorem 3.8 Let (ai(0), a2(0), • • • ) be the infinite vector of initial conditions. If the 
limit 

oo 

lim Va„(0)r" (3.2) 

r->-l+ ^—J 
n— 1 

ejc/sfs ana? is finite, then for all n € N, |a n (£)| — ► as i — > oo. 

Looking back at the examples, this result correctly separates those which decay to 
zero pointwise in n and those which do not. 

Denote by £ p the p-summable sequences: For p > 0, 

oo 

£P :={( ai ,a 2 ,...):J2K\ p <oo}. (3.3) 

n=l 

If a = (m, a2, . . •) € i 1 then (13.21 exists and is finite. For future reference, we recall 
the norm || • \\ iP defined by ||a||^ = J2™=i \ a n\ p - 

A complimentary question is to understand for which initial data the system con- 
serves energy. 

Theorem 3.9 Ifa(Q) = (a x (0), a 2 (0), • • • ) G i 2 , then 

\\a(0)\\% = ||a(i)||| 2 

for all time t > 0. 



Example 13 .21 13 . 31 and 13 .41 all have nice limits at z = — 1 in the sense of d3.2l > and 
they decay to zero as t — > oo as dictated by Theorem l3.8l It is particularly interesting 
to compare Example 13 .41 and l3~6l Theorem O . 81 correctly says that the first decays to 
zero as time increases while declining to comment on the second. 

Theorem 13.91 correctly states that Example |^] conserves energy. However, Theo- 
rem |^] is not completely satisfactory in that it only applies to solutions which have 
finite total energy. A number of our example have infinite energy. We now turn to 
understanding in detail the dynamics of dl.lt . 
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4 Solution to the initial value problem 

In this section we show that the initial value problem associated to equation II. H is 
well-posed and admits solutions for a wide class of initial data. Theorem O . 1 1 above is 
an immediate consequence of Theorem l4. 1 I below. After giving a general existence and 
uniqueness theorem, we present specific initial conditions (covered by the existence 
theorem) for which equation II. 1> admits solutions which conserve energy, dissipate 
energy, and blow up in finite time. 

We begin by formally calculating a representation of the solution given by a gener- 
ating function. We will verify that the representation is valid in the next section. Given 
initial conditions {a„(0) : a„(0) el,ne N}, we assume that a solution a n (t) exists 
and define the generating function 

oo 

G(z, t) = a n+1 (t)z n . (4.1) 

71 = 

Proceeding formally, it is straight forward to verify that G(z, t) would satisfy the fol- 
lowing partial differential equation: 

*L = tf-l)*L +jsG (4.2) 

with initial condition 

oo 

G (z) " 5>, i+1 (0)z™ . (4.3) 

n=0 

The first term on the right hand side of II. 1> comes from zj^(zG) and the second from 

One obtains an ansatz for the form of the solution by solving equation (14. 2> by the 
method of characteristics. By verifying that this ansatz solves the equation, we obtain 
the following existence and uniqueness result whose proof is postponed until section 
□ 



Theorem 4.1 Consider M.l\ with the initial condition a„(0) such that Gq(z) is ana- 
lytic in a neighborhood of the interval (a, 0]. Then the solution of M.ll exists and is 
unique for all t € [0, £*) where t* — arctanh(—a A 1). 
The unique solution is given by 

a n (t) = A ^4dz (4.4) 

with 

G(M) = ^^(GWtXz), (4.5) 
cosn(t) 

Vt(z)=- — r— - and 4> t (z) = — — . (4.6) 

1 — z tann(t) 1 — z tann(t) 

r is any simple closed contour around the origin within the region of analyticity of 
G(z, t) (which is non-empty). 



where 
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It is worth noting that G(z, t) solves the PDE given in 14.21 with G(z, 0) = Go(z) as 
initial condition. 

Notice that finite time existence of solutions only requires that the initial data a„(0) 
have at most exponential growth in n, i.e. there exists C > and 7 > 1 such that for 
all n £ N, |a n (0)| < G7™. If the a n (0) decay exponentially then the solution exists for 
all times (i.e. t* = 00). 

In addition, J4.4I implies that a n+ i(t) is the nth term of the Taylor series expansion 
of G(z, t) about z — 0. It is also straight forward to see that J4.4i defines a semigroup: 
for any suitable f(z), let 

( S tf)( z ) ■= —^r;Mz)f(Mz)) . (4.7) 
cosh(t) 

Then (I4.5i can be expressed as G(-,t) = StGo and it is easy to check that for any 

t, s > 0, 

S t ° S s f = S s o 5 t / = 5 t+s /. 

Since we are particularly interested in knowing the total energy of the solution, it is 
useful to notice that if G(x, t) is as in J4.U then 

||a(t)||| 2 J] |a n (i)| 2 = - / |G(e i9 ,t)| 2 de = ||G(e ie , t)||| 2(Sljde) . 

n=l ^ ^° 

We shall give more properties of the solutions of equation M.W in section [5] after a 
brief discussion of the examples given in the previous section. 

4.1 Analysis of examples 

We use Theorem l4.1l to calculate exact the solutions given in Section[3] 
The initial data in Example |3.2l translates in to Gq(z) = 1, so that 

G{z,t)= ' 



cosh(t) — z sinh(i) 

Calculating the Taylor series gives the quoted a n (t). In example 13.31 one obtains 
Go(z) = 1/y/l — z. Hence, 

e -t/2 1 



G(z,t) = 



Vcosh(t) V(l - z )( l ~ ztanh(t)) ' 



whose Taylor expansion produces the quoted a n (t). Example 13.41 yields Gq(z) = 
1/(1 — z), G(z,t) = e _t /(l — z ) anc l tne desired a n . Example 13.51 gives Gq(z) = 
1/Vl - z 2 and G(z, t) = - z 2 = G Q (z). Examplel3~6lgives G (z) = 1/(1 + z) 

and G(z, t) — Example l3 . 71 gives Go(z) — a/(l + az) and 

1 a 

G(z, t) = — —77- ~j j for t < t* = arctanh(l/a); 

cosn(i) — a smn(i) 1 + az 

it blows up at t = t* as stated. 
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5 Properties of the solutions 

We begin by presenting two results which are more quantitative versions of the results 
in Theorem l3.8l and Theorem l3.9l Together, they highlight the fact that it is possible to 
have a given coordinate converge to zero while no global energy dissipation is present 
in the system. This implies that there is a flux of energy out to higher and higher modes. 
It is also interesting that both of the next two results apply in some situations where the 
total energy is infinite. 

Theorem 5.1 Suppose that Gq(z) is analytic in a neighborhood of '(— 1, 0] such that 

G+(-l)= lim Go(x) (5.1) 

xeK 

exists and is finite. Then, 

a n (t) ~ 2e _t G f (-1) as t -> oo. 

In particular, in a neighborhood of the origin G(z,t) ~ 2e — *(?q (— 1)/(1 — z) as 
t — > oo. 

Remark 5.2 If a n = rT a for a < 1, then a n (f) ~ 2e _t r(l — a) as t — > oo. 

The next result contains Theorem l3.9l as well as giving control of higher Sobolev-like 
norms. We recall the Sobolev-like sequence spaces for s 6 K: 

h s = {a = (<zi,a 2 , . . .) : \\a\\ hs < oo} 

where the norm || • ||/j s is defined by ||a||? = Y^=i n2s \ a n\ 2 - 

Theorem 5.3 // 

OO 

||a(Q)||| 2 = 5^ l«n(Q)| 2 < oo 

n=l 

then 

\\a(0)\\% = \\a(t)\\% 
for all t > 0. Similarly, for any s G N, if 

oo 

||a(0)||^ 5> 2s |a„(0)| 2 < oo, 

n=l 

then for any T < oo f/zere ex/sfs a constant, C(T), such that 

sup 1100)112, < C(T). 
te[0,T] 

5.1 Finer properties of solutions 

We begin by giving conditions guaranteeing that the solution decays exponentially in 
time. We then turn to the case for which Go has singularities on the boundary of the 
unit circle. We close the section with a result which compares the dynamics obtained 
by placing a single unit of mass at different locations. 

Let D(r) denote the open disk of radius r about the origin. We begin with a simple 
criteria which guarantees that the a n decay exponentially in n. 



Properties of the solutions 



15 



Theorem 5.4 Suppose Gq(z) is analytic in the disk D(l + rj) for some r\ > 0. Then, 
for each T > 0, there exist constants 7, C > which depend only on rj and T such 
that 

sup \a n (t)\ < Cy l for all n £ N. 
te[0,T] 

In particular, the system conserves energy at all finite times. 

In order to investigate the behavior of solutions when Gq(z) is not analytic in D(l + 
rj) for some 77 > 0, we introduce the following region: 

A(C, n, 9) = {z: \z\ < ICI + V, I arg(z - - arg(C)| > 9} . (5.2) 

We begin with a careful analysis of the case when there is a single singularity on the 
unit circle. We contrast the cases for which the singularity is located at ±1 in Theorem 
I5.5l and Theorem l5.7l respectively. The remaining cases are covered by Theorem l5.10l 

Theorem 5.5 Assume that either 

a n(0) ~ Cn a ~ l as n — > 00 

for some C 7^ or, more generally, that Gq(z) satisfies both of the following condi- 
tions: 

i) There exist an A ^ and a > so that 

G a (z) - A(l - zy a fliz^l. (5.3) 

ii) Gq(z) is analytic on A(l, 77, for some rj > and < 9 < ir/2. 
Then, for all time t > 0, 

A 

a n (t) e d-2a)t n a-l ag n _^ ^ (5 4) 

F(a) 

and the solution decays to zero as t — > 00 pointwise in n; more precisely, 

a n (t) - A2 1 ~ a e~ t as t -> 00 . 

Notice that in the setting of Theorem 15.51 the energy of the system is infinite for all 
a > 0: Y^=i a n(t) = +°°- However, as a direct consequence of J5.4I > in Theorem l5.5l 
and our definition of dissipation in section|2]the following is true: 

Corollary 5.6 In the setting ofTheorem \5.5\ ifO < a < 1/2 the system is conserva- 
tive; if a = 1/2 the system displays a finite dissipation rate; and, if a > 1/2 the system 
displays an infinite dissipation rate. 

We now consider a singularity at z = — 1. The remaining points on the unit circle 
behave much like the z = 1 in that the solution decays to zero in time. They are 
discussed in Theorem 15 . 1 01 1 ater in the section. The next result shows that if there is a 
singularity at z = — 1 then the system can explode in time. 

Theorem 5.7 Assume that either 

a„(0) - C(-l)"n a_1 as n -> 00 

for some C 7^ or, more generally, that Gq(z) satisfies both of the following condi- 
tions: 
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i) There exist an A ^ and a > so that 

GoO) ~ A(z + iy a asz^-l. (5.5) 

ii) Gq(z) is analytic on A(— 1, ?y, 9)\{—l} for some rj > and < 9 < tt/2. 
Then, for all time t > 0, 

a n (t) — e (2a ' 1)t (-l) n n a - 1 as n -> oo 

r(a) 

while for n fixed 

a n (t)~2 1 ~ a e (2a - 1)t AC? l _ 1 ast^oo, 
where C% is the n th coefficient of the Taylor series expansion at z = of the function 



1 



1 



1 



1 + 2 



Notice that the energy of the system is again infinite for all a > 0: X^^Li l a r»(0| 2 = 
+oo. However, we have: 

Corollary 5.8 In the setting of Theorem 15.71 IfO < a < 1/2, f/ie system is conser- 
vative; if a = 1/2, it displays a finite explosion rate; and, if a > 1/2, it displays an 
infinite explosion rate. 



Remark 5.9 One must be careful when interpreting the results in Corollary 15 . 81 since 
the system has infinite energy. For example, when a = 1/2 then at each moment of 
time the system displays a finite explosion rate since energy is pumped in from infinity 
into any finite collection of modes. However, the rate must slow, falling to zero at 
t = oo, since as t — > oo the system converges to the fixed point given in Example l3.5l 
This follows from the fact discussed in Section ITTI that the fix point in Example 13.51 
corresponds to the initial function G(z) = Bj \/l — z 2 for some B / 0. 

We now give a more general result covering a singularity on the unit circle at any 
point other than —1. Theorem l5.5l is a special case of the following result when £ = 1. 



Theorem 5.10 Let £ be a point on the unit circle not equal to —1. Assume that Go(z) 
behaves as 



G (z) — A(( - z) a with Ay^O and a > as z — > (, 



(5.6) 



and that Gq(z) is analytic on A(C, rj, #)\{C} for some r\ > and < 9 < tt/2. Then, 
for all time t > 0, 



a n (t) 



A 



Qn a as n — > oo 



where Ct = 4>t (0- (Notice that <p t (1) = 1 and Ct — > 1 as t — > oo). The solution 
decays to zero as t — > oo pointwise in n; 



a n (t) — 2e" 



. A 
(i + O" 



as t — > oo 
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From the examples above it is natural to conjecture that any singularity on the 
unit circle dominated by a polynomial-like singularity of degree less than 1/2 will not 
destroy energy conservation. It can be shown that this intuition is correct for a wide 
class of initial conditions. We already know that if the initial conditions have finite 
energy, then energy is conserved. By the reasoning in the section defining anomalous 
dissipation, it is enough to have 

lim ni\a n \ = . (5.7) 

n — >oo 

This is possible even if the total energy is infinite. Using the Tauberian theorems in 
HHilO 1 1 we can show that the dynamics preserves a subset of sequences satisfying (15.71 . 
These solutions have infinite energy yet conserve energy in the sense of Section|2] As 
these results are tangential and a bit technical we do not give the details. 

5.2 The fixed point 

By combining the results of the previous section, one can understand a wide range 
of behavior. We illustrate this by examining the convergence to a fix point. In Exam- 
nle l3.5l we saw that the initial data corresponding to Gq(z) = 1/ VT — z 2 was invariant 
under the dynamics. This function is analytic in the open unit disk and has two square 
root singularities on the unit circle: one at z = 1 and another at z — —1. Up to a 
technicality, we show that this is the only fixed point that characterizes the initial data 
that converges to it. 

If the solution is to exist for all times the function Gq(z) must be analytic in a 
neighborhood of (—1,0]. From Theorem 15. II we see that if the limit as z — > —1 of 
Go(z) is finite then the solution converges to zero pointwise in n as t — > oo. Hence, if 
the system converges to a nontrivial steady state, it must have a singularity at z = — 1. 
We assume that the singularity is power-like (i.e. (z + 1)~ Q ). One can likely deal with 
other singularities, however, we choose not to pursue this matter here. Theorem 15.71 
implies that the singularity must be order 1/2 if it is a power; otherwise, the system 
would blow up or decay to zero. The question is: Can one make sense of the dynamics 
when \ a n\ 2 — °°? The two facts above imply that any initial condition which has 
only polynomial singularities at z — — 1 (if any) and which converges to a fixed point 
must be of the form: 

GoO) = ^-t + G (z) (5.8) 

(1+Z)2 

where A ^ and Go is analytic in a neighborhood of (—1, 0] and Go(z) — > Gq < oo 
as z — > — 1 from the right along the real axis. By Theorem l5. li the dynamics starting 
from Go converge to zero as t — > oo pointwise in the sequence space. Theorem 15. 51 
says that the first term converges to Ay / 2/(l — z 2 ) as t — > oo. Hence, all data of the 
form d5.8l > converge to the fix point from Examnle l3.5l We have proved the following 
result: 

Theorem 5.11 Let (ai(0), <i2(0), ■ ■ ■) be an initial condition such the solutions exists 
for all time and converges coordinate-wise to a fixed vector (Si, Sa, . . .) which is not 
the zero vector. Assuming that the Gq(z) associated to this initial condition has only a 
power-like singularity at z = — 1 then Gq is of the form given in J5.8i and 

°° I — 2 — 

Y j an + iz n =AJ Y — 2 . (5.9) 

n=0 

In particular, all fixed points of equation dl.lt are given by the Taylor series of (15 .91 . 
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Notice that the fixed points have infinite energy and have the property that all of the 
a n , with n odd, are equal to zero. In Section[6] we will see that steady states with less 
a degenerate structure are obtained by forcing the system. 

Remark 5.12 If one starts with an initial condition which has a polynomial singularity 
on the unit circle at £ 7^ —1 and of an order (3 + 1/2 with j3 > 0, then for all finite 
times the norm ||a(t)||/ ls will be infinite if s > —(3. Yet the dynamics still converges to 
the fix point given by Theorem l5.11l Hence, at t = 00 all || • \\h s norms are finite for 
s < 0. 

5.3 The effect of shifting the initial condition 

In Example l3.2l we described how a unit of mass placed at a\ spreads out. The follow- 
ing theorem states that the solution obtained by placing a unit of mass in a p behaves in 
the same way as the solution obtained by placing a mass at a\ except that the picture is 
shifted p units down the chain. 

Theorem 5.13 Let Gq{z) = z p for some p GN, then 

ajt) = (~l) p — h P n (t) — 

cosh(t) cosn(i) 

where £ = p — p A n and (3 n (t) satisfies 



\Pn{t)\ < 



n p [l -tanh 2p (^)] n>p 
nP tanh p ~"(f) [1 - tanh 2 "(f)] n < p 



Proof Theorem \5J3\ First observe that 



cosh _1 (Q 



G(z,t) = 



1 — z tanh(i) 
Expanding this we find 



z — tanh(i) 



1 — 2 tanh(i) 



v 



00 

G ( z > *) = H z n [(-lf tanh"+ p (i) + (3 n tanh"- p +«(i)] 

cosh(i) ' 



(3 n (t) = tanh pAn (i) [tanh-^t) - tanh(i)] k [- tanh(i)] p ~' < l ' ' < " 



n=0 

where 

pAn 

fc=i 

So 

pAn , 

< n p J2 t 1 - tanh 2 (i)] fe [tanh(i)] pAn_fc [tanh(i)] p_fe K 
fc=i ^ 

In both cases the quoted estimate follows by using the binomial theorem. □ 
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6 The forced system 

Since ( II. H may display anomalous dissipation, it is not unreasonable to expect that 
adding a forcing term to this equation may lead to a (statistical) steady state. We now 
show that this is indeed the case. Specifically, we study the system 

a n (t) = (n - l)a„_i(i) - na n+1 (t) + l n =m/(*), (6-1) 

where m € N, and f(t) is either a constant forcing term, /(f) = 1, or a white-noise 
process, fit) = W(t) (in the second case (16. 1> has to be properly interpreted as an 
infinite system of coupled Ito stochastic differential equations). 

As in the unforced setting, we represent the solution to i6.l\ as in (14.41 for some 
G(z,t). a n (t) is the nth coefficient in the Taylor series expansion of G(z,t), By 
Duhamel's principle one sees that G(z, t) must satisfy the generalization of ( 14. 21 with 
the effect of f(t)l n=m included: 

^ = (z 2 -l)^+zG + F(z,t), 
at oz 

where F(z, t) = z m ~ 1 f(t) and the initial condition is G(z, 0) = Gq{z). This equation 
is valid for both f(t) = 1 and f(t) = W(t) and forcing on any m £ N. Using the 
semigroup representation defined in J4.7I >. the solution of the equation above can be 
represented as 

G(z,t) = (S t G )(z) + [ (St-.F)(z,a)da. (6.2) 
Jo 

We have 



Theorem 6.1 Consider \6. 1\ with m = 1, f(t) = 1, and initial condition a„(0) satis- 
fying the assumptions ofTheorem \5.1\ Then, 

if n = l 
if n > 2. 



lini "»") "„ { \% r(^) (6 - 3) 



i r(^) 



Remark 6.2 Notice that 



7T 

as n 
n 



This implies that a* has infinite energy: X)^Li l a n| 2 = +°°. This is consistent with 
the fact that the steady state must be dissipative to compensate for the effect of the forc- 
ing term since dissipative solutions must have infinite energy. In this simple example 
we can see how the forcing is balanced explicitly. Mirroring the calculation in Jl .31 for 
any N: If we start the system {a* } at time t = we have 

n<N 



Hence, every block conserves energy as must happen at a fix point. 
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Theorem 6.3 Consider \6.1\ with m = 1, f(t) = W(t), and initial condition a n (—T) 
satisfying the assumptions ofTheorem \5.1\ Then, 

K - ' dW(s) a.s. (6.4) 

.oo cosh(t - s) 

Remark 6.4 From i6A\ . a**(t) is a Gaussian process with mean zero and covariance 

/ ; 2/ , {t - S) ds 

n, m € N. 



_oo cosh 2 (t - s) ( 6 « 

1 



n + m — 1 

Again, this is consistent with the need for dissipation and implies that the invariant 
measure for (16. 1> with a white-noise forcing is supported on functions with infinite 
energy (and, in particular, (16.51 is not trace-class). 

In addition, notice that this is consistent with the fact that, at least in expectation, 
the steady state needs to dissipate precisely the energy pumped into the system. In fact, 
for any N, 

\j t E Y, K* (*)| 2 = \- N®<(t)a*N + i(t) = l-^ = °- 

n<N 

Proof Theorem WT] The first term on the right hand-side accounts for the initial con- 
dition. Theorem l5. ll implies that (S^Go) — > as t — > oo. The second term is given 
explicitly by 

(S t - s F)(z, S )ds = [ A ~ S } Z) - {^(z))^ 1 f(s)ds (6.6) 
J cosh(£ - s) 



Letting f(t) = 1 and m = 1, this expression becomes 

l cosh(i — s) 



(S t - s F)(z,s)ds = I -\ s } - ds. 



It follows that 



* " ' ■ z - 1 



lim / (St- s F)(z, s)ds = ■ arctan 



a* is the nth coefficient of the Taylor series expansion at z = of this function. □ 



Proof Theorem W3\ Letting f(t) = Wit), m = 1, and considering the initial condition 
at t = — T, we have 

(S t -sF)(z, s)ds = [ ^- S } Z) dW(s). 



_ oc cosh(£ — s) 

a**(t) is the nth coefficient of the Taylor series expansion at z = of this function. 

□ 
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7 Proofs of the main theorems 

We begin by making a number of observations which will be used in the proofs. ip t 
and (f> t each have a single simple pole at z — 1/ tanh(t). Hence, at any finite time both 
are analytic in an open disk containing the closed unit disk and the Taylor coefficients 
of their expansions about zero converge to zero exponentially in n. 

For any fixed t > the map 4>t is a fractional linear transformation which bijec- 
tively maps the open unit disk onto itself and leaves the unit circle invariant. The points 
z = 1 and z = — 1 are the two fix points. In addition, for every fixed z £ D(1)\{1} and 
fixed neighborhood JV of —1, there exists a time T(z,JV) such that 4>t{z) £ Af PI -D(l) 
for all t > T(z,Af). The behavior of Go o (f> t in a neighborhood of the origin will 
be important in the analysis which follows. Observe that <j>t(0) = — tanh(i) and for 
sufficiently small r > 0, {pe tS : 9 £ [0, 2w], p £ [0, r]} is mapped approximately to 
{- tanh(i) + p(l - tanh 2 (t))e ie : 9 £ [0, 2tt], p £ [0, r]} and strictly into the closed 
disk 

E r = | - tanh(i) + p(l - tanh 2 (t))e ie : 9 £ [0, 2n],p £ [0, r/(l - r)]} . 

For sufficiently small r > 0, £ r is strictly contained in the unit disk for all times t. 
Furthermore, £ r is bounded away from the boundary by two lines emanating from —1 
of the form {-1 + pe ±w : p > 0} for some fixed 9 £ (0, tt/2). 

We recall a basic fact from complex analysis which will be used repeatedly in the 
arguments that follow. To show that the Taylor coefficients about zero of Go ° <j>t 
converges to those of F as t — > oo it is sufficient that for all t > 0, F and Go o <fi t are 
analytic in a fixed, t independent neighborhood of the origin and that Go ° <fit converges 
uniformly to F on that neighborhood. 



Proof Theorem UJA In order that equation J4.4I > be well defined, G(z,t) needs to be 
analytic in a neighborhood of z = 0. Since ip t is analytic for all z £ D(l) and each 
t > 0, the analyticity of G(z, t) about z = is equivalent to the analyticity of Go(z) 
about (f>t(0) = ~ tanh(i). As t increases, <j>t(0) decreases monotonically along the real 
axis from to —1. We need only show that Go is analytic in an open neighborhood 
of [— tanh(i), 0] in order to complete the proof that equation J4.4i is well defined for 
all t < tit. Go is analytic in a neighborhood of the closed interval [— tanh(i), 0] since 
[— tanh(t),0] C (a,0]. Appealing to the arguments stated at the beginning of this 
section we see that the image of a small ball about the origin under the mapping <j) t lies 
in a thin strip about [— tanh(i), 0]. Hence, the reconstruction formula of equation J4.4I 
is well defined because Y can be deformed to lie in a sufficiently small ball about the 
origin. 

If 1* < oo, then G(z, tanh(i*)) fails to be analytic at z = since Go(z) in not 
analytic at a = — tanh(i»); we cannot continue the solution in this case. 

To see that the a n (t) (defined as in the statement of the theorem) do define a solu- 
tion, observe that by the definition of G(z, t) and integration by parts 

. f d t G(z,f) j I (z 2 -l)d z G(z,t) + zG(z,t) j 

a n (t) = <t> dz = <J> az 

J T 2niz n J T 2mz n 

= {n-Y)j> ^iz^-i dz ~ n f 27r!z"+ 1 dZ ^^ n ~ 1 ) a «-i ^ ~ n a n +\{t)- 

This shows that a n (t) given by J4.4I > is indeed a solution of II. II for the initial condition 
a„(0) as long as G(z, t) is analytic around z = 0. 
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Finally, to show that the a n (t) defined by equation J4.4I is the unique solution 
of II. li for the initial condition a„(0), note that if two different solutions exist for 
the same initial condition, then their associated G(z, t) must both satisfy (14.2b for the 
same initial condition Gq(z). Since the solution of I4.2> is unique, this leads to a 
contradiction. □ 



Proof Theorem \5A\ Appealing to the discussion at the beginning of the section and the 
fact that cosh(f) ~ |e* as t — * oo, it is enough to show that ip t (z)(GoO(j) t )(z) converges 
uniformly to G + (— 1)/(1 — z) on some neighborhood of the origin. First, observe that 
ipt(z) converges to 1/(1 — z) as t — > oo uniformly on any disk contained within the 
unit disk. 

From the discussion at the beginning of the section we see that for all i > the disk 
of radius 1/10 is mapped to a disk contained entirely in the open unit disk and bounded 
away from the unit circle by lines emanating from —1 of a constant angle. Hence, we 
have 

lim sup \G (<p t (z))-G+(-l)\ =0. 
t ~ >DO M<i/io 

□ 

Proof Theorem \5~3\ Since (di(0), aa(0), . . .) is square summable, g(r, 9) := Go(re ) 
is in L 2 (dd) of the unit circle for all r G [0, 1]. Hence, by Plancherel's theorem, 



v />2tt -i />2rr 



|^(e t9 )| 2 
cosh 2 (/j) 



\G (<p t (e t9 ))\ 2 d9 



Introducing the change of variable e l " — 4> t (e ), we see that 



}^K(t)\ -— 2 |G (e )| — _ T —_ dry 



2^7 cosh 2 (f) |^ t Wr 1 (e*'))| 
- ^ / \G (^)\ 2 dr, = V |a„(0)| 2 . 

To obtain the bounds on the weighted norms, we notice that for s E N 

oo 

5> 2 >„(t)l 



1-2-k 




Jo 





2 

d6 . 



Using the same change of variable one easily shows that this term can be bounded in 
terms of the L 2 norms of jjgrGoie 18 ) for r < s. By assumption, these norms are finite. 

□ 



Proof Theorem \574\ Fix T > 0. By the considerations at the beginning of the section, 
one sees that for t G [0, T] the map ipt(z)(Go°4>t)(z) remains analytic in D(l+r)i) for a 
sufficiently small 771 depending only on T and -q. More precisely, 771 is picked to ensure 
that D(l + 771) maps into D(l + 77) for all t G [0, T]; the only remaining constraint on 
771 is that l + 7/i < r— t: t so that tpt(z) is also analytic in D(l + 771). Hence, the power 
series converges absolutely on D(l + 771/2) and |a„(f)| < G(l + 771/2)™ for all n G N 
for some C > 0. □ 
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Proof Theorem \5~5\ Theorem \5.1(A and Corollarv \5.6\ Theorem l5.5l is a special case of 
Theorem l5.10l so we concentrate on the later. By the discussion in the proof of Theorem 
15.11 it is clear that for each moment of time t > there exists a 771 > so that G(z, t) 
is analytic on A(£t, 771, #i)\{Ct}- "Hi mav be chosen to be sufficiently small in order to 
avoid other singularities of Go which initially lie outside A(£, 77, 6*)\{£} and approach 
D(l) under the dynamics of 4> t . A similar consideration needs to be taken into account 
for 9 and may result in an increase of 9 to a new 9%. As z — > Ct, we see that 



[cosh(i) - Ctsinh(/j)] c 



A 



[cosh(t) + Csinh(/j)] Q (Ct 
1 + C\ _ t , fl~C 



- z) a 
-1 l-2a 



A 



(Ct - zY 



G(z,t) 



The result on the asymptotics in n then follows from Theorem lA.il in the appendix 
since G(z, t) is analytic on A(Ct, T)i, Oi)\{Ct}- The result for fixed n as t — > 00 is just 
a restatement of Theorem l5.1l in this context. The Corollary follows directly from the 
discussion in section|2 □ 



Proof Theorem \5y\ and Corollarv Ulft The proof is similar to that of Theorem 15.51 
Since z = — 1 is a fixed point for <j) t for all t > 0, G(z, t) has a singularity at z = — 1 in- 
herited from Go(z). Since the circle is invariant under <j>t, for sufficiently small 771 > 
and 9\ sufficiently close to ir/2 we have that G(z, t) is analytic on A(— 1, 771, 9\). Direct 
calculation yields: 



G(z,t) 



1 



cosh(7j) + sinh(i) 



1 + tanh(t) 
1 - tanh(t) 



A 



(1 + z) a (1 + z) c 



as 



-1 . 



We obtain the quoted result by applying Theorem lA. 1 I from the appendix. 
The asymptotics in time follow from the fact that 



G(z,t) 



Ae (2a-l)t / X 



1 



1 



as t 



00 



and direct expression of the right hand side in a Taylor series in z. Corollarv l5 . 8l follows 
from the discussion on anomalous dissipation in section|2]and the above results. □ 



8 Inviscid limits 

We return to the analysis of the inviscid limits of II. 11 1. Fixing peN and defining 

p 

A„ = Y[(n k) 

k=\ 

with the convention that A n = 1 if p = 0, we consider 

d„,„ = ~2vA n a nil/ + [(n - l)a„_i^ - na n+ i, v ] + l n= iW(t) . (8.1) 

As mentioned in Section fTTTl it is straightforward to see that this system converges to a 
random variable a** = (a\* v , a^, • • • )■ In fact, one has 

E^A„|<J 2 = i. 

n 
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Thus, the system does not display anomalous dissipation; the dissipation which bal- 
ances the energy injection (due to the forcing) comes from the term — 2vri p a n (t). 
Setting 

OO 



one sees that 



= (z 2 - 1)^ + zQ v - 2uz"^- + W{t) . (8.2) 
ot oz az p 

Using the variation of constants formula we obtain 

Gui;t) = S+ v g v i-,G)+ f (5 t _ s „1) dW s . 



ii 



We will concentrate on the case p 6 {0, 1}. By the method of characteristics, we find 
that 2 ^ 

\^77-^t(z)f(Mz)) forp = 

iS t , v f)iz)=\ cos \® 

7-—^^ K ,t(z-v)f(v + K(f> K ,t(z-v)) forp=l 

I coshint) 

where k 2 = 1 + v 2 . 

It is interesting to contrast the regularizing effect of the different terms. When 
p = 0, St.v simply dissipates energy at a faster rate than St- When p = 1, St, v has 
a stronger regularizing effect than St^, in that the characteristics are attracted to the 
circle (y — n)e lB inside of the unit disk and the singularity of ip K tiz — v) stays uniformly 
bounded outside of the unit disk for all times. Hence if / has a radius of convergence 
greater than v — k ^ — 1 + v then S"' 1 f in analytic on a disk with radius greater than 
one all times uniformly. 

For fixed t, St^f converges to St as v — ► uniformly on a neighborhood of the 
origin. Since one also has that Stf, St^f all go to zero uniformly on the open disk as 
t — ► oo for / bounded on the unit disk, we have that J_ (<St_ S)1/ )l dW s converges to 

f (St- a l) dW s in mean squared as v — » 0. As before we are primarily interested in 
these solutions. In this setting they are given by: 



cosh(i — s) 

k tanh(K(i — s)) 



(83) 



-dWis) p = 1 



j [k + v tanh(K(£ — s))] n cosh(K(i — s)) 
Theorem 8.1 For alln GN, t € E, and p — 0, 1, one has 

lim E[<*„(t) - <*(t)] 2 = 

u— >0 1 

Hence a** v (f) converges to a**(i) as v ^ for p = 0, 1. Furthermore, one has the 
estimates given in (II . 14-i one/ dl.l5t . 

Proof Theorem WT\ . Fix any i. Consider a„(i, i>, 0) and a n ii) starting from initial con- 
dition zero at time T with T < t. As T — > — oo, a n it, v, 0) and a„(i) converge 
respectively to a**(i, v, 0) and a**(i). 
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By the same argument as Theorem l6.3l one see that J8.3I > holds. Subtracting J8.3i 
from (16. 4>. one obtains 



~,2(n-l) 

\a* n *(t, v, 0) - a;;(t)} 2 = 2 I [(1 + uT v - l \ \ u + 2) 2n du 



< 



r°° i 



(1 + + 2i/)' 



which implies that a**(i, i/, 0) — * a**(i) almost surely as i/ — * 0. The convergence of 
in the other cases is similar. Applying the Ito isometry to 18.31 proves the quoted value 
of E[a**(t, v, p)} 2 forp = 0. The other estimates follow from 

1 1 1 

< < 



k + v k + v tanh(t) k 
which holds for t > 0. □ 

Remark 8.2 At first glance, it might seem more natural to consider the system 

a n ,v = -2vn p a n<ll + [(n - l)a n _i )V - na n+ {\ + \ n =iW(t) . 
This leads to the following equation for the generating function Q v : 

^ = (z 2 - + zQ v 2vz^ - 2UT>» + W(t) . 

at oz oz p 

where V p is p applications of the operator defined by (Vf)(z) — §^(zf(z)) and V° is 
the identity operator. Hence, we see that the extra dissipative term contains derivative 
of all orders less then or equal to p. Not surprisingly, the result is a mixture of the 
behavior of ( 18. II for all orders less than or equal to p. In particular, when p = 1 the 
asymptotic (in time) behavior is given by 

/"* Ke~ v{t - S) tanh( K (i - s))"- 1 

- J ^ ^ K + v tanh ( K ( t _ s ))jn cosh(K;(t - s)) (S) 

and satisfying the following estimate: 

9 A second linear shell model 



We begin the analysis of the second model 11.1 li by giving general conditions for the 
existence of a unique solution of the initial value problem. The technique is similar to 
that used in Section^] 

Theorem 9.1 Let {b n (0)} be such that 

oo 

^(-l) n 6 n (Q)< oo. (9.1) 
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Then the solution of HI. Ill exists and is unique for all positive times. It can be repre- 
sented as 

b " (t)= (2~n-l)! ft^-i ^Q (9 ' 2) 
where t 

H(x, t) = E x H (X(t)) exp ( - 1 J X 2 (s)ds) , (9.3) 

anof X(t) satisfies the stochastic differential equation 



dX(t) = -X 3 (t)dt + -^=y/l- X\t)dW{t). 
v2 



(9.4) 



E x denotes the expectation conditional on X(0) = x G [—1,1] and 

H (x) = ^(-ir +1 6„(0)x 2 "- 1 . (9.5) 

Remark 9.2 Alternatively, H(x, t) can be expressed as 

H(x,t) = ^(-lr+V^z 2 "- 1 . (9.6) 

riGN 

where 6„(i) solves ( 11.1 H . 

Remark 9.3 If the sequence {b n (0)} is monotone and converges to zero as n — > oo 
then the condition in ( I9.lt holds. 

The following theorem summarizes the most interesting properties of solutions 
of (fTTTl. 



Theorem 9.4 Suppose that b n (0) satisfies i9.lt . Then, for any positive time t > 0, 

lim (2n + l)& 2n +i(*) = C\(t) and lim 2nb 2n (t) = C 2 (t), (9.7) 

n — >oo n — >oo 

where C\(t),C2(t) € K, Ci(t), C^it) ^ Q for all but finitely many t S [0, oo). /n 
particular, there exists a T > smc/i f/iaf /or aZZ t > T, the solution of il.llt is 
dissipative and satisfies 

OO OO 

n— 1 n— 1 

/n /acf, ^ZnLi &n(*) — * as t — > oo. 

Remark 9.5 The fact that equation ( II. lit dissipates energy at finite times is implicit 
in the representation ( I9.3I >. As time grows, the factor exp(— i J Q X 2 (s)ds) converges 
to zero as exp(— ct) almost surely for some positive deterministic c. (This follows from 
the law of large numbers and the verifiable assumption that the process is ergodic.) 
Hence, H(x, t) converges to zero uniformly in x as t — > oo. 
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Writing ( 19.31 as H(x, t) — {T t Ho)(x), it is easy to see that T t defines a (Feller) 
semigroup with generator L defined by 



18 (<'-**>£)-?'■ 



(9.9) 



4<9x 

for / £ C 2 ([-l, 1]). In addition, H(x, t) satisfies 



dH _ 1 d f„ ^dH\ x 2 
~dt ~ 4dx 



with initial condition H(x,0) — Hq(x) for x £ [—1, 1]. One can check that the 
boundaries at x = ±1 are entrance boundaries for (19. lOt and H(x, t) satisfies 

lim (1 - x 4 )— = 0. (9.11) 
x->±i ax 

Proof Theorem WT\ . Noting that 

Lx 2 "- 1 = -n(n + \)x 2n+1 + (n - l)(n - ±)a; 2 "~ 3 , 

we compute 

oo 

^(-D" +1 6„(i)a ; 2 ^ 1 

n=l 

oo 

= ^(-l)' i+1 ((n - l)(n - i)6 n -iW - n(n + §)6»+i(t)) x 2 "" 1 

n=l 

oo 

= ^(-l) ,l+1 6„(<) (-n(n + i).T 2,i+1 + (n - l)(n - |)z 2n - 3 ) 

n=l 
oo 

= ^(-l)"+ 1 6 n (f)L(ar 2n - 1 ). 



n=l 



□ 

Proof Theorem WM . Associated with J9. lQi we have the eigenvalue problem 




(9.12) 



subject to the boundary conditions 

lim (1 - x 4 )^- = 0. 
a:-»±i ax 

It is straightforward to see that the operator in (19.121 equipped with the boundary con- 
dition in ( 19. lit is self-adjoint in L 2 [— 1, 1]. We now explain why this operator has 
discrete spectrum. A standard calculation shows that the boundary is an "entrance 
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boundary" in the sense of Feller (|Fel54 McK56|), i.e. the diffusion ( I9.4> . if started 
from the boundary, enters (—1,1) and does not return to the boundary. 

Define L v = \d x (l - x 4 )^ + V(x) where V(x) = -^.By standard PDE theory, 

dtu — L v u, 

subject to the condition lim^o u(t, x) = S y (x), has smooth solution in (—1,1) for 
any t > 0. We denote this solution by pY(x,y). For a fixed t > 0, p((x,y) is 
Lipschitz for x E [—1 + e, 1 — e] and y S [—1, 1] with a fixed Lipschitz constant C t e 
and sup.,, », e [_i i] \pY( x > V)\ < D t . Consider the solution to the following initial value 
problem: Let / G L 2 ([— 1, 1]) and solve 

dtu = L v u with it(0, x) = f(x). 

The solution is given by 

u(t, x) is Lipschitz for x S [— 1 + e, 1 — e] as the following simple estimate shows. 

\u(t, x) — u(t, x')\ = 

< Ct\x - x'\ J \f(y)\dy < 2CI\x - x'\ [J \f(y)\ 2 dy 
u(t, x) is also bounded in terms of ||/||2 as follows: 



l 

(pY( x , V) - Pt(x', y))f(y)dy 
i 



i 

• / \pY(x,y)- P Y(x',y)\\f(y)\dy 



\u(t, x)\ 



J ^ pY(x, y)f(y)dy < D t J ^ \f(y)\dy <2D t (J^ |/(y)| 



As one can see by a Cantor diagonalization argument in intervals I n = [— 1 + 1 — 

i], T t = e tL is a compact self-adjoint operator. Therefore, the spectrum of L v is 
discrete. 

Note that the lowest eigenvalue has the following variational representation: 

r Shi - ^'(x)? + \x 2 cj> 2 (x))dx 
A = mi j , 

<t> S-i <t> 2 (x)dx 

where the infimum is taken over L 2 equipped with the boundary conditions ( 19. 111 . 
This shows that the spectrum is strictly positive. Let {<fik(x), Xk}ken be the pair of 
eigenfunction and eigenvalues such that each 4>k(x) is odd in x (the even ones do not 
matter since the initial condition Hq(x) of ( 19. lOt is odd from ( 19. 5i ). The solution 
of ( 19 . 1 Oi can be represented as 

H(x,t) = J2 h ke~ Xkt 4>k(x), (9.13) 

where 

l 

H (x)4>k(x)dx. 
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In turn, ( 19. 131 implies that 



where p\ is defined by 



b n (t) = J2 h ^ Xkt Pn 



Aes 



(9.14) 



(9.15) 



The p„ satisfy the following recurrence relation inherited from (11.1 II : 

-\ kP k n = (n - l)(n - - n(n + \)p k n+l , n e N, p£ - 0. 

The following lemma describes the asymptotic behavior of np n . 



Lemma 9.6 for every > 0, f/ze recurrence relation in \9.16\ implies that 
lim (2n + V)p\ n+1 = c\ and lim (2n)p 2n = ci 

where c\ and cf are nonzero constant whose sign is the same as that ofp\. 
Proof. Assume p\ > and write ( 19. 161 as 



Pn+l 



/ i 1 \Pn 

n(n + y) 



n(n + 5 ) 



(9.16) 



(9.17) 



We omit the index k in this proof as it plays no role. For sufficiently large n and C 
depending only on A, 



Pn+l < 1 



2 C 



max{p„,p„_i}. 



This implies {p n } is bounded. More is true: 



Pm < 



J—n 



max{p„,p„_i}. 



This implies that 



mp m < exp|logm+ i ^ log ^1 - y + } max{p n ,p„_i j 

which implies further that lim sup mp m < oo. On the other hand, 

2 



Pn+l > 1 - 



n + 1 



Pn-1, 



which implies 



z + i 
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This implies that 



mp m > exp | log to - i ^ log ^1 - j-^jj | Pn if n-m=0 mod 2. 

Thus, Urn inf mp m > 0. To show that the the sequences in the theorem are Cauchy 
simply compute 



\(n + l)p„+i - (71 - l)p n -l\ = 



A 



n(n + 1) 



(n + l)p n+ i +(n- l)p n . 



1 1 



2 n(n + ±) 



Using the fact that {np n } is bounded in n and summing over n we see that the sequence 
is Cauchy and have proved the lemma. □ 



Going back to the proof of Theorem l9.4l using (19. 17i in J9.14I implies (I9.7I > with 



Ci(t) = e~ Xkt h k cl and C 2 (t) = £ e'^^c; 



keN 



Note that there is a T > so that Ci(t)C 2 (f) > for all t > T. Finally, implies 
that 

1 d °° 

xtE = - J im N < N + §)M*)&tf+i(*) = -^i(t)^ 2 (t) < 

z at ^ — ' too 



for all i > T which proves 

We also consider the system of forced equations 

b n = (n- l)(n - |)6„_i - n(n + hb n+ i + f(t)l n=m , 



□ 



(9.18) 



for n = 1,2,... with boundary condition &o(t) = for all t and / (t) is either a constant 
forcing term, fit) — 1, or a white-noise process, fit) = Wit). 
We have 



Theorem 9.7 Consider A9. 181 with fit) — 1, and initial condition 6„(0) satisfying 19.il ). 
77ien 

t — >oc 



lim Mf) = b* n = V 
where pj? is defined by \9.15\ and 



A-ei 



cf>j s: (x)x rn dx 

In particular, 6* satisfies 

lim (2n + 1)&2„ +1 = CJ > and lim (2n)6^„ = C 2 * > 0. 



Theorem 9.8 Consider \9.18\ with fit) = Wit), and initial condition b n (—T) satis- 
fying i9.ll . Then 



lim b n (t) = b* n *it) = V d k p k n f 

fcGN J ~ 



t- XkS dWis) 



where p„ is defined by \9.15\ 
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In particular, b**(t) is a Gaussian process with mean zero and covariance 

E6«(*)0*)= E ^tfy^ 

fe,fc'GN fe 

and we have 

lim (2n + l) 2 E(&o« + 1(*)) 2 = C** > and lim (2n) 2 E(6^(t)) 2 = C 2 ** > 0. 

n— *oo n— *-oo 

Appendix A Estimates on Taylor coefficients 

For the reader's convince, we now state a theorem on the asymptotic of Taylor's series 
which can be found in |FO90). . 

Theorem A.l Let A(£, 77, 6*) fee 05 in (15.21 . Assume that f(z) is analytic in A(£, 77, #)\{C} 
/or some ( e C, i) > 0, ant/ < < 7r/2. T/' 

f(z) ~ as z — > C 

/or some if > anc/ a ^ {0, — 1, —2, • • • } f/ien 



T(a) c™ +a 

where f n is the n-th Taylor coefficient of f(z) about z = 0. 
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